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ABSTRACT 

Based on high-resolution ultraviolet spectroscopy obtained with the Far Ultraviolet Spectroscopic 
Explorer (FUSE) and the Cosmic Origins Spectrograph, we present new detections of O vi and N v 
emission from the black-hole X-ray binary (XRB) system LMC X-3. We also update the ephemeris 
of the XRB using recent radial velocity measurements obtained with the echelle spectrograph on the 
Magellan-Clay telescope. We observe significant velocity variability of the UV emission, and we find 
that the O VI and N V emission velocities follow the optical velocity curve of the XRB. Moreover, the 
O vi and N V intensities regularly decrease between binary phase = 0.5 and 1.0, which suggests that 
the source of the UV emission is increasingly occulted as the B star in the XRB moves from superior to 
inferior conjunction. These trends suggest that illumination of the B-star atmosphere by the intense 
X-ray emission from the accreting black hole creates a hot spot on one side of the B star, and this hot 
spot is the origin of the O vi and N v emission. However, the velocity semiamplitude of the ultraviolet 
emission, K\jy rs 180 km s^" 1 , is lower than the optical semiamplitude; this difference could be due 
to rotation of the B star. Comparison of the FUSE observations taken in November 2001 and April 
2004 shows a significant change in the O vi emission characteristics: in the 2001 data, the O vi region 
shows both broad and narrow emission features, while in 2004 only the narrow O vi emission is clearly 
present. Rossi X-ray Timing Explorer data show that the XRB was in a high/soft state in the Nov. 
2001 epoch but was in a transitional state in April 2004, so the shape of the X-ray spectrum might 
change the properties of the region illuminated on the B star and thus change the broad vs. narrow 
characteristics of the UV emission. If our hypothesis about the origin of the highly ionized emission is 
correct, then careful analysis of the emission occultation could, in principle, constrain the inclination 
of the XRB and the mass of the black hole. 
Subject headings: stars: individual (LMC X-3) 
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1. INTRODUCTION 

How do black holes affect and interact with their sur- 
roundings? This is one of the highest priority questions 
of current astrophysics. Typically, the angular momen- 
tum of accreting material drives the formation of a disk 
around the black hole. The hot inner accretion disk is 
a powerful source of high-energy X-ray radiation, since 
immense gravity can cause substantial heating of matter 
spiraling in toward the black hole. In general, the X-rays 
emerging from the inner disk and corona can illuminate 

1 Based on observations with (1) the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
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sion, operated by Johns Hopkins University, supported by NASA 
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scopes located at Las Campanas Observatory, Chile. 
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and ionize surrounding material, and in the case of black 
holes in binary star systems, the high-energy emission 
can affect matter on the surface of the companion star, 
in the stellar wind, in the accretion flows, or the outer 
disk. A significant amount of the X-ray radiation could 
be reprocessed into ultraviolet/optical photons. Such re- 
processing might account for the continuum and/or some 
of the observed emission features in the UV/optical spec- 
tra of black-hole systems (although UV/optical contin- 
uum emission can also emanate directly from the accre- 
tion disk). For this reason, the study of UV spectra of 
black hole systems can provide valuable insights about 
the physics of accreting black holes and the impact they 
have on their surroundings. The analysis of the variabil- 
ity of the UV radiation from a binary star system as a 
function of orbital phase can determine the location of 
the region producing the UV photons, the physical con- 
ditions existing in the region, the mechanism producing 
the radiation, and kinematic parameters of the binary 
itself. 

LMC X-3, which was first discovered by ILeong et al.l 
(J1971I) , is an excellent target for this type of study. It is 
a bright X-ray binary (XRB) system in the Large Mag- 
ellanic Cloud composed of a B3 V star (V ~ 16.7 — 17.5) 
and a central black-hole (Cowley et al. 1983, hereafter 
C83). It is one of the few black- hole candidates that are 
persistently luminous in X-rays and also bright enough 
to be observed with high-resolution ultraviolet spectro- 
graphs with good signal-to- noise (S/N). The black hole is 
believed to be undergoing accretion from its B-star com- 



Song et al. 



panion via Roche lobe overflow with an orbital period of 
1.7 days (C83; Soriaet al. 2001). Spectroscopic observa- 
tions of the B star indicate a large radial velocity semi- 
amplitude, K B = 235±11 km s' 1 (C83). Although LMC 
X-3 has been intensively studied in the X-ray band, UV 
spectroscopic observations of this target are more sparse. 
The first Far Ultraviolet Spectroscopic Explorer (FUSE) 
observation of LMC X-3 was made by Hutchings et al. 
(2003, hereafter H03) in November 2001. During the 24 
ks exposure, LMC X-3 was in its brightest X-ray phase. 
These observations revealed a broad emission feature in 
the vicinity of the O vi AA1031.9, 1037.6 doublet, which 
Hutchings et al. argue is due to the blend of the 1031.9 
and 1037.6 A lines (see their Figure 4). From this emis- 
sion, they measured a shift in the O vi emission velocity 
between two binary phases of about 100-150 km s^ 1 . As- 
suming the O vi emission arises in the inner parts of the 
black hole's accretion disk, they conclude that the mini- 
mum mass of the B star and the black hole are 13 and 15 
Mq, respectively. However, their measurement of the ve- 
locity variation of O vi emission with binary phase suffers 
from some uncertainties. First, the ephemeris they used 
to determine the phase is from more than 20 years ago. 
Over such a long time baseline, the errors contained in 
the ephemeris could be accumulated, resulting in phases 
that are off by «0.1 at the la level. Another source of 
uncertainty is the limited phase coverage of their obser- 
vations. Due to the modest S/N of their FUSE data, 
Hutchings et al. summed their FUSE spectra into only 
two phase bins centered at phases = 0.53 and 0.70, which 
only loosely constrains K^h, the velocity semiamplitude 
of the black hole. Overall, these limitations preclude any 
firm conclusions about the nature and implications of the 
O vi emission. 

An updated ephemeris and expanded UV spectroscopy 
with better orbital coverage are needed to better con- 
strain the implications of the UV emission lines and the 
characteristics of this XRB system. In this paper, we 
revisit LMC X-3 with a new ephemeris and new UV ob- 
servations. We update the ephemeris of LMC X-3 based 
on high-resolution optical spectroscopy recently obtained 
with the 6.5m Magellan-Clay telescope. We also report 
new ultraviolet observations of LMC X-3 from two high- 
resolution instruments. First, we present temporal moni- 
toring of the O VI emission using FUSE data with a much 
longer time baseline that provides more than three quar- 
ters of the LMC X-3 orbit. With the extended FUSE 
data, we detect narrow O vi emission from the binary 
system and measure the variations in velocity and in- 
tensity of the emission as a function of orbital phase. 
Second, we complement the O vi analysis with new ob- 
servations of the N v AA1238. 8, 1242.8 doublet obtained 
with the Cosmic Origins Spectrograph (COS) on board 
the Hubble Space Telescope (HST). The N v emission is 
detected at high significance in the COS data and pro- 
vides corroborating evidence of the velocity and inten- 
sity variations of the highly ionized emission. During the 
time of the FUSE observations in 2004, LMC X-3 was 
also observed quasi-simultaneously in the X-ray band- 
pass with Chandra and the Ross i X-ra y Timing Exper- 
iment (RXTE); see iWang et all (|2005l ) for full details. 
While our emphasis in this paper is on the UV emission 
lines, we also briefly comment on the affiliated X-ray ob- 



servations. With the combination of this information, 
we investigate the following questions: (1) Does the O vi 
emission originate in the accretion disk, the ionized stel- 
lar wind, the illuminated surface of the B star, or a hot 
spot (e.g., from the spot where overflowing Roche lobe 
material impacts the accretion disk)? (2) Does the N v 
doublet emission arise in the same region as the O vi 
emission, i.e., are the N v and O vi features physically 
associated? 

2. OBSERVATIONS 

2.1. Optical Spectroscopy 

Fifty three spectra of the optical counterpart (Warren 
& Penfold 1975) of the X-ray binary system were ob- 
tained on the nights of 2005 January 20-24, 2007 Decem- 
ber 20-21, and 2008 February 27 - 2008 March 1 using 
the M agellan Inamori Kyoce ra Echelle (MIKE) spectro- 
graph (Bernste in et al.ll2003l ) on the 6.5 m Magellan-Clay 
telescope at Las Campanas Observatory. The instru- 
ment was used in the standard dual-beam mode with a 
l'.'O x 5'.'0 slit and the 2x2 binning mode. With a few ex- 
ceptions, the typical exposure time was 1800 s. We only 
focus on spectra from the blue arm, which has a wave- 
length coverage of 3430-5140 A and resolving power of 
R = 33,000. The radial velocities were determined from 
these new spectra using the FXCOR task within IRAF 
and then fitted to a circular orbit model. The typical 
uncertainties of the radial velocities are 3.5 km s . 

2.2. Ultraviolet Spectroscopy 
2.2.1. FUSE Observations 



Foll owing the initial observations of iHutchings et al.l 
(12003[ ). LMC X-3 was observed again with FUSE in April 
2004, and during the latter campaign quasi-simultaneous 
observations were recorded with Chandra and RXTE 
(jWanget al.ll2005f ). The 2004 FUSE observations pro- 
vided a total exposure time of 86 ks («4x longer than 
the original FUSE observations obtained in H03). LMC 
X-3 is a variable source, and in April 2004 the average 
UV flux of the target turned out to be « 4 x fainter in 
the r elevant spectral r e gions compared to the UV flux 
when IHutchings et al.l (|2003t ) observed it in November 
2001. Consequently, the S/N ratios are comparable in 
the FUSE data obtained in 2001 and 2004. However, the 
addition of the later data set significantly improves the 
orbital coverage of the O Vi data. Both the 2001 and 
200 4 FUSE data were reduced with CALFUSE, version 
3.0 (|Dixon et al.ll2007l ). It is important to note that the 
FUSE SiC channels suffer from scattered solar light prob- 
lems in the 2004 data, and in this paper we only use the 
data from the LiF channels (details regarding the design 
and p erformance of FUSE can be found in IMoos et al.1 
l2lMl200l . 

In Figure [TJ we compare the FUSE observations of 
LMC X-3 obtained in 2001 (black histogram) and in 2004 
(red histogram) in the wavelength range near the O vi 
doublet (in this figure, all of the data obtained in 2001 
are accumulated into a single spectrum, and likewise all 
of the 2004 data are coadded). Bright airglow emission 
lines from the Earth's atmosphere are marked (©), and 
well-detected absorption lines from the O VI doublet and 
C II A1036.3 near v (heliocentric) « km s _1 are also 
labeled (analysis of the absorption lines in the FUSE 
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Fig, 1. — C omparison of the FUSE observations of LMC X-3 obtained in 2001 (black histogram, see H03) and in 2004 (red histogram, see 
Wang ct al. 2005) in the vicinity of the O VI A1031.9.1037.6 doublet. In this figure, all data from the November 2001 observing campaign 
have been accumulated into a single spectrum regardless of the phase of the binary during the individual exposures, and likewise all of the 
April 2004 observations have been coadded to form a single spectrum. The O VI doublet emission, whose velocities and intensities vary 
with the binary phase as presented in Figures [4] and \E\ are smeared out in the combined spectra. Both spectra are binned to 7.5 km s 
pixels. Airglow emission lines from the Earth's outer atmosphere are marked (©), and absorption lines due to O VI and C II (A1036.3) are 
also labeled. The weaker O VI absorption line at A1037.6 is blended with Galactic absorption due to C II* (A1037.0) and H2 (A1037.2). It 
might also be blended with C II at the velocity of the LMC itself. 

this region can only be observed with COS at low spec- 
tral resolution {R sa 2000). COS could observe LMC 
X-3 with higher resolution (R w 20, 000) at A > 1130 A, 
but at the expense of missing the OVI doublet. How- 
ever, other highly ionized species are accessible in the 
A > 1130 A range (e.g., N v AA1238.8, 1242.8 or C iv 
AA1548.2, 1550.8). 

We have recently observed LMC X-3 with COS using 
the high-resolution G130M and G160M modes (Program 
11642). The target was observed on 15 December 2009 
with total exposure times of 16.9 ks (G130M) and 2.7 
ks (G160M). The data were calibrated with the pipeline 
CALCOS (version 2.11b), and major fixed-pattern noise 
features due to grid wires in the cross delay line detector 
were removed with a flat field developed by K. France. 
This flat field does not correct fo r all types of fixed - 
pattern noise in COS spectra (e.g.. lSahnow et al~ll2010H . 
but for the measurement of broad and bright emission 
lines presented in this paper, the lower-level fixed pattern 
noise does not have a significant impact on the results. 

The regions of the COS G130M and G160M spectra of 
LMC X-3 covering the N v and C iv emission lines are 
shown in Figure [21 Emission in t hese species w a s pre - 
viously detected from LMC X-3 bv lCowlev et all (|1994f ) 
using low-resolution spectra obtained with the Faint Ob- 
ject Spectrograph on HST as well as the International 
Ultraviolet Explorer. Our new COS data offer two ad- 
vantages over the previous observations: (1) the spectral 
resolution is «10 times higher, so variations in the ve- 
locity centroids can be measured more accurately, and 
(2) the S/N is higher, and consequently the data can be 
divided into phase bins for investigation of emission vari- 
ability vs. binary phase. As we will show below, the N V 



spectrum of LMC X-3 is reported in I Wang et"alll2005h . 

The substantial variation in the UV flux is immediately 
obvious, but we note that the flux variation cannot be 
described with a simple scale factor applied to the entire 
spectrum. For example, the continuum is roughly two 
times brighter at A « 1045 A but is roughly 4x brighter 
near the putative O VI emission at A « 1034 A. It is dif- 
ficult to draw conclusions about the long-term UV flux 
variability with observations at only two epochs, but it 
appears that while the O VI emission intensity decreased 
along with the continuum in the 2004 observation, there 
is not a strict correlation. We will show below that nar- 
row O VI emission lines are detected in the 2004 data, but 
it is possible that the broad O VI emission was not present 
in the later observation for some reason. The shape of 
the X-ray spectrum of LMC X-3 also changed between 
2001 and 2004 (see H3.4[) . so it is possible that the changes 
of the broad O VI feature are related to changes in the 
X-ray flux. 

2.2.2. COS Observations 

There are some limitations of the spectroscopic UV 
dataset from FUSE presented in the previous section. 
First, LMC X-3 is a challenging target for FUSE because 
the UV flux is relatively low, so the data can be noisy. 
Also, in the FUSE wavelength range, emission lines from 
the Earth's atmosphere are a potential source of confu- 
sion (see below). Moreover, FUSE is no longer operating, 
so we cannot contemplate follow-up FUSE spectroscopy. 
Fortu nately, all of these problems can be overcome with 
COS ([Green et al.H2003D . the new spectrograph on HST, 
with some caveats. COS has limited sensitivity in the 
912 < A < 1130 A range (jMcCandliss et all 1201(1 . but 
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Fig. 2.— Portion of the COS G130M (left) and G160M (right) spectrum of LMC X-3 (black histogram) showing the emission lines of 
the N V A1238.8.1242.8 doublet and the C IV A 1548.2,1550.8 doublet. The N V and C IV doublets are both apparent in emission as well 
as absorption. Although these features seem reminiscent of P-Cygni profiles, it is more likely that the absorption features are due to the 
absorption by the foreground ISM and the emission lines are from LMC X-3. See Section 3.3 for more detail. In the left panel, the smooth 
black curves show the model of the underlying LMC X-3 flux used to extract the N V emission and absorption characteristics. The red 
dotted curves show only the fitted N V emission lines, and the green curves show the total N V fits (absorption and emission). The COS 
data are binned to 10 km s — 1 pixels in this figure. 

riod. The periodogram for the combined data is shown 
in the middle panel of Figure [3] The minimum value 
of xl occurs for P = 1.7048089 ± 0.0000011 days and 
T = HJD 2, 454, 454.9964 ± 0.0011, which are very close 
to the values derived from the MIKE data alone (top 
panel of Figure [3]) but with much smaller uncertainties 
in P. In this paper, we use this more accurate ephemeris 
to determine the orbital phase of the UV observations of 
LMC X-3. The bottom panel in Fig. [3] is the zoom-in 
near the minimum \t of the middle panel. This panel 
clearly shows that all other possible alias periods are 
ruled out with high confidence. The fitted velocity am- 
plitude is 250.3 ± 1.1 km s _1 , comparable to the value 
derived by Cowley et al. (1983). A detailed descrip- 
tion of the optical spectroscopy and a full discussion of 
the improved dynamical model used here to obtain the 
radial velocities will be presented in a separate paper. 

3.2. O VI Doublet Emission 

The FUSE observations cover about three quarters of 
the XRB orbit. In order to measure the O vi emission 
velocities and intensities as a function of orbital phase, 
we divide the FUSE exposures into 5 binary phase bins, 
which provides reasonably good resolution of trends vs. 
phase while still maintaining adequate S/N in the spectra 
from each bin. Fig. |4] shows the spectra of these 5 phase 
bins. 

In this and several other figures, we refer to the 2001 
FUSE observations from H03 as the "old" data and the 
later FUSE observations as the "new" data. The old data 
were binned in orbital phase in a similar way and are 
also represented in Figure |U Due to the low-Earth orbit 
of the FUSE satellite, many airglow emission lines from 
the Earth's outer atmosphere are pres ent in the FUSE 
wavelength range (Fel dman et al.ll2001[) . Strong airglow 
emission lines are readily apparent and easily identified 
(marked with ©), but weaker telluric emission could be 
a source of contamination. Most of the airglow emission 
lines are excited by sunlight and are only present in data 
recorded on the day side of the FUSE orbit. Since the 
photons are time tagged, it is straightforward to extract 
spectra using only photons detected during the night por- 



velocities are consistent with the optical orbital veloci- 
ties of LMC X-3 and are also consistent with the veloc- 
ity trends exhibited by the O vi doublet. The G160M 
spectrum is considerably noisier and the C iv emission 
features are weaker, so in this paper we will focus on the 
velocity and amplitude variations of the N V lines and 
their comparison to the O VI emission. 

3. RESULTS 
3.1. New Ephemeris of LMC X-3 

We have analyzed the optical observations of LMC X-3 
in a manner analo gous to the recent study of LMC X-l 
(jOrosz et al.l 120091 ). Among the 53 optical spectra that 
we obtained with the Magellan-Clay telescope, we ex- 
clude seven observations from our analysis due to their 
relatively large residuals in the radial velocity^ The re- 
maining 46 radial velocities are then fitted to a circu- 
lar orbit model, which returns the systemic velocity Vo, 
the time of inferior conjunction of the optical star Tq 
(phase 0), and the velocity semiamplitude of the sec- 
ondary Kb, when given the value of the orbital period 
P. We construct a periodogram as shown in the top 
panel of Figure [3] using a wide range of trial periods. 
The best fit indicates P = 1.704820 ±0.000012 days and 
T = HJD 2, 454, 454.9962±0.0011, where the uncertain- 
ties on the individual measurements are scaled to give 
1. 

We also tried including the radial velocities published 
in C83 in order to refine our measurement. However, 
the systemic velocities obtained from the two data sets 
are quite different: 406 km s _1 for MIKE and 310 km 
s _1 for Cowley et al. The cause for this difference is 
not clear, although it is not uncommon to find that 
the systemic velocity changes depending on what line 
is used (one even sees shifts when one uses different 
Balmer lines). We thus first simply subtracted the re- 
spective systemic velocities from the two data sets be- 
fore x 2 fitting sine curves as a function of trial pe- 

9 Five of the seven outliers were obtained on the same night, so 
it is possible that a subtle problem with the wavelength calibration 
on that particular night compromised those measurements. 
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Fig. 3. — Top: The periodogram of LMC X-3 derived from the 
Magellan+MIKE radial velocity measurements, constructed by fit- 
ting a three-parameter sinusoid at each trial period. The y-axis 
shows the value of xt of the fits- Middle: The periodogram de- 
rived from the combination of the MIKE measurements and the 
Cowley et al. (1983) radial velocities. Bottom: An expanded view 
of the middle panel near the minimum xi ■ 



tion of the orbit. While this reduces the overall S/N, 
night-only spectra are useful for supressing the strong 
airglow emission lines (which sometimes blend with - 
or mask entirely - spectral features of interest) and for 
assessing whether an emission feature is a bogus identifi- 
cation due to confusion with terrestrial airglow emission. 
So, we also show in Figure|4]some spectra extracted using 
night-only data. Since the velocities of the emission lines 
are of particular interest, we also show in Figure [5] the 
phase-binned FUSE spectra plotted versus heliocentric 
velocity. 

The spectra in Figures 0] and [5] are arranged from top 
to bottom according to phase (as indicated by the labels 
in each panel). From inspection of these figures, we can 
recognize emission from both lines of the O vi doublet 
that regularly shifts in velocity as a function of the XRB 
orbital phase. Interestingly, the intensity of the O vi 
emission also appears to change with phase: the O vi 
emission appears to be brightest near inferior conjunction 
of the X-ray source (orbital phase « 0.5), and it fades and 
disappears into the noise at superior conjunction (phase 
w 0). As we will show below, similar behavior is observed 
in the N v doublet. To quantify this observation, we have 
fitted Gaussians to the emission features, and the O vi 
velocity centroids and line widths indicated by the fitting 
exercise are summarized in Table [1] 

As noted above, one concern about the identification 
of the O VI emission is that the features could be con- 
taminated by terrestrial airglow lines - there are sev- 
eral molecular N2 emission features in the vicinity of 
the O Vi A1031 . 9 feat ure, e.g., N 2 b(l,2) and N 2 b'(9,6) 
(jFeldman et alj 120011 ). and strong telluric emission is 
present near O vi A1037.6. However, there are several in- 
dications that these emission features we identify as O VI 
emission lines do not arise in the Earth's atmosphere. 
First, the velocity centroids of the features regularly shift 
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(a) (new) all data: phase 0.37-0.55 




(d) (new) all data: phase 0.56-0.70 




(e) (new) night only: phase 0.56-0.70 

(f) (old) all data: phase 0.65-0.7 





(i) (new) all data: phase 0.0-0.20 
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Fig. 4. — Orbital modulation of the O VI emission lines. Each 
panel shows the spectrum obtained by coadding all observations 
in the phase range indicated in the panel (e.g., panel a shows the 
data accumulated when the XRB was between orbital phase = 
0.37 and 0.55). Panels labeled as "new" data are based on obser- 
vations obtained in 2004, and panels labeled as "old" show data 
from the 2001 observations. "All data" indicates that the spectrum 
includes data obtained during both the day and night sides of the 
FUSE orbit and thus is affected by terrestrial dayglow emission 
lines (marked with ©). In some cases the O VI A1037.6 transition 
is significantly blended with terrestrial emission so we also show the 
spectrum obtained by combining only the nighttime data (which 
greatly suppresses the airglow emission features). In some panels, 
Caussian fits to the O VI emission is overplotted with a smooth red 
line. As discussed in the text, the O VI emission lines show regular 
variations in velocity and intensity as a function of orbital phase. 



with orbital phase; this would not occur if the lines were 
telluric. Second, there are stronge r N 2 features at othe r 
wavelengths in the F USEba,nd (see lFeldman et al"1 [2001). 
and the absence of those stronger N 2 features indicates 
that the N 2 b(l,2) and N 2 b'(9,6) emissions are negligi- 
ble. Third, in the night-only spectra the airglow emis- 
sions in this range disappear but the candidate LMC X-3 
emission lines remain apparent. Finally, the O vi trends 
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TABLE 1 
LMC X-3 O vi Emission Line Measurements 



O vi 

Transition(A) 



FUSE 
Observation Date 



Binary 
Phase 



Vo vi 

(km s~ x ) 



FWHM 
(km s~ 1 ) 



1031.93 


2004-04 (day+night) 


0.37-0.55 


385 ± 11 


171 ± 26 


1031.93 


2001-11 (day+night) 


0.42-0.52 


384 ± 13 


126 ± 31 


1037.62 


2004-04 (night only) 


0.37-0.55 


350 ± 13 


143 ± 32 


1031.93 


2004-04 (day+night) 


0.56-0.70 


218 ± 9 


106 ± 22 


1037.62 


2004-04 (night only) 


0.56-0.70 


207 ± 8 


59 ± 16 


1031.93 


2001-11 (day+night) 


0.65-0.71 


154 ± 8 


53 ± 16 


1031.93 


2004-04 (day+night) 


0.70-0.83 


147 ± 8 


68 ± 18 
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Fig. 5. — Same data as shown in Figure HI but plotted vs. helio- 
centric velocity. The zero velocities are centered at Galactic O VI 
absorption lines (1031.93 A and 1037.62 A), respectively. The cen- 
troids of the emission features associated with LMC X-3, as deter- 
mined from Gaussian fits to the lines, are indicated with red tick 
marks. 



are corroborated by the N v emission lines in the COS 
data, which are clearly detected at high significance and 
are not affected by telluric lines. 

3.3. N V Doublet Emission 

The region of the COS G130M spectrum of LMC X- 
3 shown in Figure [2] clearly shows significant emission 
from both lines of N v along with weaker N v absorption. 
Upon initial inspection, the N v features seem reminis- 
cent of P-Cygni profiles (the C iv doublet shows similar 
P-Cygni-like characteristics), and it is tempting to at- 
tribute the features to a stellar wind from the B-star 
companion of the LMC X-3 black hole. While this is a 
possible origin for the emission features, one must bear in 
mind that the foreground highly ionized absorption from 
the ISM of the Milky Way (and from the corona of the 
LMC ) is strong in the direction of the Mage llanic Clouds 
(e.g.. (Wakker et al J 11991 IHowk et all 12003) . so the N v 
and other highly ionized absorption is likely to arise (at 
least partially) in the foreground ISM and might have 
nothing to do with LMC X-3 itself. We also note that 
the LMC X-3 spectral features appear to be narrower 
than the P-Cygni profiles typically seen from normal B 
stars (cf. iGradv et al.lll987t iPrinia 1 119891 ). and the vari- 
ability of the emission suggests an alternative origin, as 
we discuss below. 

The COS observation of the LMC X-3 covers the XRB 
orbit between phase 0.66 and 0.89 based on our new 
ephcmeris. In a fashion analogous to our treatment of 
the FUSE data, we divide the G130M exposures into 
five phase bins and plot the velocity variation of both 
lines of the N v doublet with binary phase in Figure [6j 
We used Gaussian and Voigt functions to model, respec- 
tively, the N v emission associated with LMC X-3 and 
the N v absorption arising in the foreground Milky Way 
gas. In each phase bin, the centroids and widths of the 
Galactic N v absorption lines were determined by using 
all the observations (instead of only the subset data for 
that particular single phase bin) since any foreground 
Milky Way absorption should not change with the bi- 
nary phase. As can be seen from Figure [6l this proce- 
dure produced reasonably good overall fits, which sup- 
ports the notion that the absorption is mostly due to 
foreground material unrelated to LMC X-3. It is impor- 
tant to note that while the core of the COS line-spread 
function (LSF) provides the exp ected spectral resolution , 
the COS LSF has broad wings (|Ghavamian et alJ feOQg). 
To account for these wings, the model profiles were con- 
volve d with the COS L S Fs, at appropriate wavelengths, 
from Ghavam ian et al.l (|2009D . The velocity centroids 
and line widths of the N v lines measured in this way are 
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Fig. 6. — Velocity variation of the N V emission with phase, 
plotted vs. heliocentric velocity. The zero velocities are centered 
at the rest-frame wavelengths of the doublet (the 1238.82 A and 
1242.80 A lines are shown in the left and right panels, respectively). 
The thick black lines indicate the effective continuum placement 
adopted for the N V fits, the red dotted curves show the emission 
line profiles only, and the green curves show the overall fit including 
N V absorption and emission. The spectra were rebinned to 10 
km s _1 pixels for this figure. 

summarized in Tabled 

As we can see from Figure [6] and Table [21 the velocity 
centroids of the N v emission lines do not change much 
in the five phase bins extracted from the COS data set. 
However, the COS observations were randomly sched- 
uled with respect to the XRB phase (i.e., no scheduling 
constraints were placed on the COS program in order 
to observe the system at some particular phase), and as 
we will show in §??, the COS observations happened to 
be obtained near an extremum of the XRB radial veloc- 
ity curve when the velocities are not expected to change 
much. Moreover, the N v velocity centroids are con- 
sistent with the O vi velocities recorded in roughly the 
same phase (see below). Perhaps more interestingly, the 
N v intensities change with phase in the same way as the 
O vi: the N v emission intensity decreases steadily going 
from phase « 0.65 to phase rs 0.9 (see Figure |6|). This 
is surely an important clue about the origin of the N v 
and O VI emission, and we will return to this in §??. 

3.4. X-ray Observations 

Figure [JJ compares the X-ray characteristics of LMC 
X-3 during the FUSE observations in 2001 and 2004. Al- 
though LMC X-3 usually occupies a soft state in which 
the X-ray spectrum is dominated by emission from the 
accretion disk, it was in a faint and transitional state 
near the time of the 2004 FUSE observations (this is in- 
dicated by the presence of a s trong power-law component 
in the observed spectrum; see lCui eta l. 2002). As shown 
in Figure El LMC X-3 appeared to be nearly as bright in 
X-rays during the 2004 FUSE observations as it was dur- 
ing the 2001 observations of H03. This is not accidental; 
the 2004 FUSE observations were part of a coordinated 



Chandra/FUSE target-of-opportunity program that was 
triggered when the RXTE ASM showed that LMC X-3 
was entering a bright state[J3 However, as can be seen 
from the left panel of Figure [71 the 2004 increase in X- 
ray brightness was associated with a short-lived X-ray 
flare, whereas the 2001 data were recorded when LMC 
X-3 was undergoing one of its more gradual (and more 
typical) upward modulations in X-ray brightness. Thus, 
while the X-ray brightness was comparable in 2001 and 
2004, it is likely that the object was in a different type of 
X-ray state. To further compare the X-ray characteris- 
tics of the XRB during the FUSE observations, the right 
panel of Figure [7J shows spectra from pointed RXTE ob- 
servations obtained very close to the times of the FUSE 
studies. From this comparison we see that LMC X-3 was 
in the usual soft state around the time of the 2001 obser- 
vations and the associated spectrum is well fitted with 
a multicolor disk model. However, LMC X-3 was in a 
significantly harder spectral state when the 2004 FUSE 
data were obtained. The spectrum for the latter obser- 
vation requires a second power-law component aside the 
disk component. The best-fit parameters from the two 
models are summarized in Table [31 

4. DISCUSSION: NATURE OF THE UV LINE EMISSION 

Where do the ultraviolet emission lines originate? 
Emission from O vi arising in interstellar plasm a has 
been frequently detected (e.g., IDixon et' al. 2006), but 
strong N v emission lines like those shown in Figure [2] 
are never seen from interstellar clouds, so it is most likely 
these emission lines are intrinsic to the LMC X-3 system. 
Moreover, lines from interstellar gas should be station- 
ary, but the LMC X-3 emission-line velocities seem to 
follow the velocity curve of the XRB. 

To show the UV emission-line velocity trends more 
quantitatively, we plot the velocity centroids of the O vi 
and N v emission lines as a function of the binary phase 
(derived from the new ephemeris) in Figure [51 In this 
figure, the O vi and N v emission velocities are plotted 
with red and blue filled circles, respectively; the horizon- 
tal bars on each point indicate the phase range associated 
with each measurement, and the vertical bars represent 
the formal ±1<7 velocity uncertainties. When we fit a 
sinusoidal curve to the UV data (green curve), we find 
an orbital semiamplitude of 180 ± 1 km s _1 and a phase 
deviation (compared to the optical data) of 0.027±0.004. 
The systemic velocity of LMC X-3 inferred from the fit 
to the UV measurements is 301 km s _1 , comparable to 
the systemic velocity (310 ± 7 km s _1 ) obtained in C83, 
assuming a circular orbit. The cyan dashed-dotted line 
shows a fit to the UV velocities in which the sinusoid 
is forced to be in phase with the optical curve but the 
semiamplitude is allowed to vary; we see that this fit is 
comparably good. Therefore, the UV emission velocities 
generally follow the trend of the optical velocities in the 
orbital phase, but the UV velocity semiamplitude is lower 
than the optical semiamplitude (250.3 ± 1.1 km s _1 ). 

If the O vi and N v emission lines arise in the atmo- 
sphere of the B star, they should follow the optical orbital 
velocities of the star. For a circular orbit, which is most 

10 This was done to maximize the S/N of the Chandra grating 
spectroscopy, which was originally initi ated f or the detection of 
X-ray absorption lines. Sec Wane ct al. (2005) for further details. 
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TABLE 2 
LMC X-3 N v Emission Line Measurements 



N V COS Observation Binary V N v FWHM 

Transition(A) Start Time(UT) Duration(s) Phase (km s _1 ) (km s" 1 ) 



1238.82 


2009-12-15 04:07:56 


1885 


0.66-0.68 


142 ±4 


96 ±7 


1242.80 


2009-12-15 04:07:56 


1885 


0.66-0.68 


140 ±3 


101 ±8 


1238.82 


2009-12-15 04:47:49 


3014 


0.68-0.70 


134 ±4 


94 ±7 


1242.80 


2009-12-15 04:47:49 


3014 


0.68-0.70 


131 ±3 


113 ±8 


1238.82 


2009-12-15 07:57:18 


5411 


0.76-0.80 


125 ±4 


87 ±10 


1242.80 


2009-12-15 07:57:18 


5411 


0.76-0.80 


126 ±3 


45 ±5 


1238.82 


2009-12-15 10:40:19 


3270 


0.82-0.85 


145 ±6 


70 ±11 


1242.80 


2009-12-15 10:40:19 


3270 


0.82-0.85 


140 ±5 


44 ±8 


1238.82 


2009-12-15 12:22:19 


3270 


0.87-0.89 


136 ±11 


113 ±23 


1242.80 


2009-12-15 12:22:19 


3270 


0.87-0.89 


124 ±10 


46 ±9 
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Fig. 7.— Left: X-ray light curve of LMC X-3 from the RXTE All-Sky Monitor with the dates of the "old" (2001) and "new" (2004) FUSE 
observations marked (dashed lines). The two insets show the regions in the vicinity of the two epochs with a time span of 350 days. Right: 
RXTE Proportional Counter Array X-ray spectral energy distribution of LMC X-3 obtained near the times of the two FUSE observations. 
The solid lines are the best-fit models to the data (see Table [3] for model details). The grey dashed and dotted lines indicate the disk and 
power-law components used to fit the new observation, respectively. 



TABLE 3 
Best-fit Parameters for RXTE Spectra of LMC X-3 a 



Observation 
Time 


(10 20 cm" 2 ) 


Disk Component 6 
kT dbb 

(keV) 7V dbb 


Power-law 

r 


Component 


2001 
2004 


4.74 
4.74 


-, 9C -+0.01 
1,zo -0.01 

1 iQ+0.01 
1 - ly -o.oi 


27.85±g : p 

24.45+^ 


2 30+ 015 
z - ou -0.16 


0.085i° ; o i 9 8 



a The uncertainties shown represent 90% confidence intervals. 
b Tdbb an( l -^dbb are th e temperature at inner disk radius and the normalization 
(dimensionless) in the multi-black body model (diskbb in XSPEC), respectively. 
c r and N po are the photon index of power law and the normalization (photons 



keV cm ' 
respectively. 



s at 1 keV) in the power law model (powerlaw in XSPEC), 



likely for LMC X-3 since the 1.7 d orbit should quickly 
circularize, the radial velocity curve should be a stan- 
dard sinusoid. This behavior is approximately seen in 
the UV emission velocities shown in Figure [8] It is pos- 
sible that there are systematic uncertainties in the UV 
velocities (e.g., due to line shapes that are more complex 
than our assumed single Gaussian plus interstellar ab- 
sorption) that are not adequately reflected by the formal 
error bars from the single-Gaussian fits, but the velocities 



are sufficiently well-constrained to show that they follow 
the general velocity trend of the optical B-star measure- 
ments and that the UV semiamplitude is clearly different 
from the optical semiamplitude. 

The similar (but not identical) behavior of the UV vs. 
optical velocities could naturally occur if one side of the 
B star is illuminated and hence heated by the X-ray/UV 
continuum emission originating near the black hole. In 
this case, the UV and B-star orbital velocities would gen- 
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Fig. 8. — Velocity variation with binary phase for the O VI (red 
filled circles) and N V (blue filled circles) emission lines. No mea- 
surements are available for the O VI emission lines in the phase 
range 0.83-0.20 since the emission is too weak to be detected near 
the X-ray source superior conjunction. Horizontal and vertical bars 
on each point indicate the phase coverage and ±1<t velocity uncer- 
tainties associated with each measurement, respectively. A sinu- 
soidal fit to the UV data is presented as the green curve. The cyan 
dashed-dotted line shows a fit to the UV data in which the sinusoid 
is forced to be in phase with the optical curve but the semiampli- 
tude is allowed to vary. The systematic velocity of LMC X-3 from 
the fit to the UV data (301 km s" 1 ) is marked as the dotted black 
line. 



E 
o 

en 

CD 



i 
o 



CD 
O 
i_ 
cn 
CD 





4 
1 


\ 






12 


VI 1031.93 (new data) • - 
VI 1037.62 (new data)* - 
VI 1031.93 (old data) 


8 


- 




- 


4 




h 


* 


n 






- 




- 


*♦ 






3 




N V 1238.82* - 
N V 1242.80* - 






~ 


2 


- 


*# 


- 








• 


1 


- 




• :* j 


n 






• 



0.4 



0.6 



0.8 



1.0 



orbital phase 



Fig. 9. — Variability of the integrated flux of O VI and N V 
emission as a function of orbital phase. 



its orbit, and since the velocity centroid would be some 
sort of weighted mean of the portion of the spot visible 
to the observer, the effective component of the rotation 
velocity along the line of sight could be less than the full 
v sin i value. To test this hypothesis, it would be helpful 
to obtain observations of LMC X-3 covering the phase = 
0.2 — 0.5 portion of the orbit. It would also be helpful 
to carry out detailed modeling of this scenario, but such 
modeling is beyond the scope of this paper. 

However, we note that the intensity variations of the 
O vi and N v emission lines are broadly consistent with 
this idea. Figure [9] shows the variations of the integrated 
flux in the O vi and N v emission vs. orbital phase. 
The gradual decrease of the integrated O vi and N v 
flux occurs over binary phases when the accreting black 
hole moves from inferior conjunction to quadrature to 
superior conjunction. At superior conjunction, the O vi 
emission is too weak to be reliably measured, and the 
N v emission is greatly suppressed (see Figures [5] and 
[6]). During this time period, the hemisphere of the B 
star directly facing the black hole would move from fully 
illuminated to partially covered to completely occulted. 
The observed UV emission intensity variations appear to 
be consistent with this hypothesis. 

A similar hypothesis has been proposed to explain the 
similar N v emission lines and other emission features 
detected in t he UV spectra of a low-mass X-ray binary 
Cygnus X-2 (jVrtilek et al.l 12003). By measuring the ra- 
dial velocities of the emission lines and comparing the 
observed line profiles with the predictions from models of 
line emission from an X-ray-heated accretion disk corona, 
Vrtilck et al. (2003) suggest that most of the line emis- 
sion detected in the UV spectra of Cygnus X-2 is from 
the illuminated surface of the companion star. For LMC 
X-3, if the idea is supported by future work, this could 
provide an important constraint on the inclination of this 
X-ray binary (and thus the mass of the black hole) as well 
as insight on how black holes affect their surroundings, 
and the properties of their companions in particular. 

Our interpretation of the detected variable O VI and 
N v emission from the X-ray binary LMC X-3 seems 
promising, but it is far from being conclusive. There 
exists considerable uncertainties in decomposition of the 
N v emission feature (§3.3). There could be weak P- 
Cygni profiles in the N V and C IV regions. For the ab- 
sorption part of the P-Cygni profile, the importance and 
shape of the wind contribution versus the contribution 
from the foreground ISM is difficult to distinguish and 
accurately measure. We also can not fully exclude con- 
tributions from other plausible components of the LMC 
X-3 system, such as stellar and/or accret ion disk winds 
as detected or predicte d in other systems (Boroso n et al.1 
l2007t iRaymond 1119931 ). Our hypothesis requires further 
observational and theoretical testing. 



erally follow each other, but rotation of the B star could 
cause the UV emission to have a somewhat different semi- 
amplitude. Cowley et al. (C83) report that the B-star 
has v sin i — 130 ± 20 km s _1 . A hot spot on the B 
star that remains pointing in the direction of the black 
hole would only be rotating toward or away from the ob- 
server depending on where the star is in the course of 
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